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1 INTRODUCTION 



nary and multiple systems (e.g. 


Duauennov & Mavoil 199ll: 
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Dal I2OO5I; 
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Kouwenhoven et al. 2005 


, I2OO7I: lLadall2006: Zinnecker & Yorkd 


20071; iGoodwin et al.ll2007r). It is also thought that the majority of 



ABSTRACT 

Over the past few decades, numerous wide (> 10^ au) binaries in the Galactic field and halo 
have been discovered. Their existence cannot be explained by the process of star formation or 
by dynamical interactions in the field, and their origin has long been a mystery. We explain the 
origin of these wide binaries by formation during the dissolution phase of young star clusters: 
an initially unbound pair of stars may form a binary when their distance in phase-space is 
small. Using A^-body simulations, we find that the resulting wide binary fraction in the semi- 
major axis range 10^ au < a < 0.1 pc for individual clusters is 1 - 30%, depending on the 
initial conditions. The existence of numerous wide binaries in the field is consistent with ob- 
servational evidence that most clusters start out with a large degree of substructure. The wide 
binary fraction decreases strongly with increasing cluster mass, and the semi-major axis of the 
newly formed binaries is determined by the initial cluster size. The resulting eccentricity dis- 
tribution is thermal, and the mass ratio distribution is consistent with gravitationally-focused 
random pairing. As a large fraction of the stars form in primordial binaries, we predict that 
a large number of the observed "wide binaries" are in fact triple or quadruple systems. By 
integrating over the initial cluster mass distribution, we predict a binary fraction of a few per 
cent in the semi-major axis range 10^ au < fl < 0.1 pc in the Galactic field, which is smaller 
than the observed wide binary fraction. However, this discrepancy may be solved when we 
consider a broad range of cluster morphologies. 

Key words: Binaries: general - star clusters - methods: A^-body simulations 



It is well known that binaries are dynamically processed in 
star clusters with wider and less bound systems t ending to be de- 
stroyed by encounters ( lHeggielll975l : lHillslll975l) . Therefore, the 
field binary population is dynamically processed wi th respect t o 
the birth population of binaries (Rroupa 1995; Parke r et alj|2009l) . 
The origin of most field binaries can be understood as a mixture of 
differently processed initial populations jGoodwinll2009l) . 



majority of binarie|3 in the field population presumably originate 
from clustered star formation. 
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' For brevity we will use 'binaries' to mean 'multiples' of any multiplic- 
ity for the remainder of this paper, only drawing a distinction where it is 
necessary. 



However, a significant number of very wide (a > 10^ au) bi- 
naries have been observed in the field (see § O. As such wide bi- 
naries are extremely sensitive to destruction they have been used 
to constrain the properties of the Milky Way. Wide binaries in the 
Galactic disc and halo have been used to place limits on the density 
of MA CHOs and other unseen material (e.g,. 'Bahcal l et all 19851 ; 
lOuinn et al. 2009), to constrain the formation history of the Galaxy 
(e.g., Allen et al. 2007), and to test the dark matter hypothesis (e.g., 
iHemandez & Leea2008,) . But the extreme sensitivity to destruction 
makes the survival (and even formation) of such binaries in a clus- 
ter something of a mystery. 
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Many very wide binaries iiave separations comparable to the 
average interstellar separation in clusters (typically a few 10'' au), 
and the very widest binaries have separations of order the size of a 
young cluster core (typically a few 10"* au). Given this, it is difficult 
to s ee how they could even form, let alone survive, in a cluster (see, 
e,g.. lScallvetai]|l999l J Parker et al ]|200^ . Even in an isolated star 
forming region the typical siz e of a star forming core is only 10* au 
dWard-Thompson et alj20 0'?) which presumably sets the very max- 
imum size of a (primordial) binary system. 

It is possible that a wide binary forms via dynamical interac- 
tions in the Galactic field (the capture mechanism). A prerequisite 
for this mechanism to work is that a significant amount of kinetic 
energy is dissipated. This energy dissipation can occur due to tidal 
friction and due to three-body interaction. Tidal friction occurs in 
the rare event of two stars nearly colliding in a close encounter. In 
the vast majority of the cases this results in a merger or a fly-by, 
and only in a small number of cases does this lead to the formation 
of a binary system. However, all binaries resulting from capture by 
tidal friction are very tight, with orbital period of several days. 

Another possible mechanism is three-body interactions. In this 
case the third star acts as the energy sink, and is generally ejected 
with high velocity. However, the stellar density in the field is low, of 
order 0.1 Mq pc"'', so that three-body enc ounters are rare, and cap- 
ture by dynamical friction rarely occurs. iGoodman&Hu^ jl993l) 
find that the creation rate Ng for binaries per unit volume can be 
approximated by 



Nb = 0.75 



(1) 



where M is the typical mass of a star in the field, n the number den- 
sity of stars, cr,, the velocity dispersion, and G the gravitational con- 
stant. In the solar neighbourhood M a 0.3 Mg, n » 0.03pc"' and 
cr,, a 50 kms"'. For the field therefore, Ng 4x 10"^' pc"' Gyr"'. 
This shows that the formation of binaries in the field is extremely 
rare. Note that in dense star clusters the stellar density is high 
and the velocity dispersion modest, such that the number of bi- 
naries formed via three-body interactions (Eq. [T]l may be substan- 
tial. However, wide binaries, which have semi-major axes compa- 
rable to the size of these clusters, are not formed, as they simply 
do not fit in these sta r cluste rs. Furthermore, W-body simulations 
by iKroupa & Burker3 ( 1200 ll) have shown that the observed broad 
period distribution of binaries in the field cannot be produced by 
dynamically modifying a tighter period distribution in a star clus- 
ter. 

A third possibility for the origin of wide binaries is formation 
during cluster dissolution, which is the mechanism we propose in 
this paper. In an evolving star cluster, stars that are initially un- 
bound, may become bound to each other as the cluster expands, 
i.e., if the gravitational influence of the other cluster members de- 
crease^ In order to form a binary pair in this way, (i) the two stars 
need to be sufficiently close together, (ii) the two stars need to have 
a sufficiently small velocity difference, and (iii) the newly formed 
binary should not be destroyed by gravitational interaction with the 
remaining cluster stars or field stars. 

Throughout this paper we refer to binaries with 10' au < a < 
O.I pc as the wide binary population. The paper is organised as 



^ Interestingly, Levison et al. (2009) have independently proposed that 
large populations of comets may be captured by stars during cluster dis- 
solution, by a mechanism whic h is similar to that di scussed in this paper for 
wide binary formation (see also lEggers et al.ll997l) . 
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Figure 1. The observed semi-major axis distribution, f(a), of the 
wide binary populat i on (in arbitrary units). The curves indicate the 
iDuquennov & MavoJ jl99lh dist ribution (assuming a mass of 1 Mp for 
each b inary), and the results from'Close et al. 1 1990), Lenine & Boneiorncl 
i2007h (for both the Galactic disk and halo). .Chaname & Gould t2004.) . and 
iPoveda et al] )2007l) . 
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Figure 2. The cumulative distributions in period an d semi-major axis for 
solar-t ype stars in the solar neighbourhood, for the iDuquennov & MayoJ 
)l99lh log-normal period distribution. The total mass of each binary is as- 
sumed to be Mr = 1 Mg. The dotted lines indicate the median semi-major 
axis, and a = lO' au. 



follows. In §[2] we briefly discuss surveys of wide binary systems 
and the corresponding observational techniques. In §[3] we explain 
our technique and assumptions. In §|4]we provide analytical and 
Monte Carlo estimates for the resulting wide binary population, 
and in § [5] we present estimates based on W-body simulations of 
evolving star clusters. Finally, in § |6] we present and discuss our 
conclusions. 



2 OBSERVATIONS OF WIDE BINARIES 

Observations have indicated that binaries as wide as I pc exist in the 
halo, while in the Galactic disc the widest binaries have separations 
of order 0.1 pc (e.g.. Close et al. 1990; Chaname & Gould 20o4). 
Wider binaries are rare, although some authors claim evidence for 
binary and higher - order mu ltiple syst ems wider than O. I pc (e.g., 
ISchoIz et aI.ll2008l:ICabalier o 2009; M amaiek et ai]|2009l) . 

Statistical properties of a wide binary population are of- 
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ten recovered using the ang ular two-point correlation function 
(e.g. , Bahcall & Soneira 1981;i Garnaviclj|l988l : lGould et al.lll995l : 
iLon ghitano & Binggeli 2010). The most prominent disadvantages 
of this method are the inability to identity individual wide bina- 
ries and the need for a good model for the stellar population that is 
studied. 

Individual wide binary candidates are often ide ntified by their 
comm o n proper motion on t he sky ( e.g., Wasserma n & Weinberg 
' 3 ll 



199ll: IC haname & Gould |2004| - lL6pine & Bongiornd 12001 : 



Makarov e t al. 2008, and numerous others). Many w ide binaries 



were f ound by Hipparcos i IeSAI 19971 ) as well; see also Soderhjelm 
( |2007|) . Their nature can then be further c onstrained by mea- 
suring the parallax and radial velocity (e.g., iLatham et aLlll984l : 
IHartigan et al.lll994lOuinn et alj2009f) . 

The most well-known survey for binarity is arguably that of 
iDuguermov & Mayor! ( 1 199 ih . who carried out a large binarity study 
for solar-type stars in the solar neighbourhood. They found a log- 
normal period distribution f{P) with a mean (log/') = 4.8 and a 
standard deviation o-\a^p = 2.3, where P is the orbital period in 
days, in the range 1 < P < 10'" days. Note that the latter value 
roughly corresponds to an orbital period of 30 Myr. In this log- 
normal period distribution, ~ 15% of the binaries have a semi- 
major axis wider than 10'' au (see Fig.O. 

Several observational studies suggest a semi-major axis distri- 
bution of the form f{a) oc which corresponds to a flat di stri- 
butio n in logg, also known as Opik's law (Opik 1924; van Albada 



1968| ; IVereshchagine t al. 1988; All en et al.ll200d ; |Poveda & Alle 



2004h . In particular, iPoveda et aL (i2007h find binaries in the field 



follow Opik's law in the separation range 100 < a ^ 3000 au, and 
suggest that a population of very young binaries follows Opik's law 
up to a K 45 000 au (0.2 pc). 

Many other authors have also found significant wide bi 



nary populations, notably Close et al. l ll99d) , lLepine & Bongiomol 
( 12007 ). and C haname & Goul d ( 20040. We summarise the wide bi- 
nary separation distribution from various authors in Fig.[T] 

The reliability of the observed properties of wide binary pop- 
ulations remains an issue. It is difficult to confidently establish 
whether the two components of a candidate wide binary are truly 
bound, or whether it is merely a chance superposition. Due to the 
long orbital periods of wide binaries, up to millions of years, it 
is impossible to accurately derive orbital properties and therefore 
confirm the bound state of the candidate wide binary. It should be 
noted that due to the lack of detailed orbital information, the quoted 
separations of wide binaries are usually the instantaneous angular 
separation. The true separations may thus be significantly larger 
than the observed, projected separations. On the other hand, an es- 
timate for the semi-major axis distribution of an ensemble of binary 
syste ms can statistically be obtained from the projected separation 
(e.g.. lLeinertetalJl993l) . 

On the other hand, it is extremely difficult to identify binaries 
with distant and faint stellar or substellar compani ons, due to con- 
fusion with foreground and background stars (e.g., IChandrasekliaJ 
[l944). The wide binary fraction as identified in the observational 
papers may therefore be a lower limit, rather than an upper limit. 

Although the exact form of the semi-major axis distribution 
for very wide binaries is not yet known (see Fig. [TJ, two features 
appear to be clear: (i) the binary fraction in the separation range 
10'' au < a < 0.1 pc is roughly 15%, and (ii) a sharp drop-off in 
the separation range 0.1 - 0.2 pc is present, likely due to dynamical 
destruction of the most weakly bound binary systems. 
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Table 1. Properties of the models used in this paper. The quantity R de- 
scribes the virial radii of models PI and P2, while it describes the radius of 
the sphere that includes the fractal structure for models Fl and F2. 



2.1 Stability of very wide binaries 

Whether a wide binary in the Galactic field is stable or not, 
depends primari ly on its semi - major axis a. Using a Monte 
Carlo approach, IWeinberg et al] jl987h show that binaries with 
a = 0.1 pc are able to survive in the Galactic disk for ~ 10 Gyr, 
roughly the age of the Galaxy, but they do not find a sharp 
cut-off in the semi-major axis distribution. On the other hand, 
a sharp drop-off is observed for binaries at a a 0.2 pc due to 
interactions with other stars, molecular clouds, and the Galac- 
tic tid al field (e.g,. Bahcall & Soneira 1981 ; iRetterer & Kind 
1982'; Mallada & Fernande^ 120011 : IJiang & Tremaind 120091 ). 
Gould & Eastman (2006) explain the slope-change at ~ 3000 au 



in the separation distribution of iLepine & Bongiomol ( |2007|) as 
the result of dynamical interactions in the Galactic field. Very 
low-mass binary systems are more weakly bound than their 
solar- mass analogues; their typical separation is therefore expected 
to be considerably smaller (e.g., iBurgasser et al.u2003i , see also 
iKraus & Hillenbrand l2009bl) , although several very low mass 
binaries with separa tions up to ~ 0.1 pc have been detected (e.g., 
iRadigan et al.ll2009l and references therein). 



3 METHOD AND ASSUMPTIONS 

Given the apparent difficulty of forming extremely wide binaries in 
star clusters, and the even greater difficulty in keeping them bound 
we suggest that these extremely wide binaries form during the dis- 
solution of a cluster into the field. 

Man y young star clusters do not survive for more than a 
few Myr (ILada & Laddl2003l ; iFall et al.|[200^ ; iMengel et al.ll2005l ; 
iBastian et alj|200^ . Their rapid destmction is probably due to the 
expulsion of the residual gas left-over after s tar formation, which 
dram a tically changes the cluster potential dGoodwin & Bastiai] 
l200d ; lGoodwidl2009l . and references therein). Stars that are un- 
bound in a cluster potential may become bound to each other af- 
ter dissolution as the local density decreases, thus forming an ex- 
tremely wide binary. 

In this section we describe the cluster models used in the re- 
mainder of this paper. We use two approaches to investigate the 
formation of wide binaries during cluster dissolution: Monte Carlo 
simulations (§|4j2} and /V-body simulations (§|5](. In § |3.1| we ex- 
plain our choices for the models used in our analysis, in § 13.21 we 
define the quantities we use to describe binarity, and in § 13.31 we 
describe the algorithms we use to ensure the stability of the newly 
formed wide binaries. 
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3.1 Model setup 

We simulate star cluster s using the STARLAB package 
dPortegies Zwart et alJ 1200 ih . We draw A' single stars from 
the iKroupj j200lh mass function, /m(M), in the mass range 
0.1 < M < 50 Mq. The lower-limit corresponds to the hydrogen- 
burning limit. The upper limit is (somewhat arbitrarily) set to 

50 Mq. 

We perform simulations with varying A', ranging from small 
stellar systems (or sub-clumps) with A' = 10 to open cluster-sized 
systems {N = 1000). We additionally perform simulations of clus- 
ters with different radii (0.1-1 pc), to identify the relation between 
the initial cluster size and the properties of the newly formed wide 
binary population. We study two sets of dynamical models: Plum- 
mer models and substructured (fractal) models, which we describe 
below. The properties for the subsets of models are listed in Table[T] 

(i) Pliimmer models. The Plummer sphere is often used in star 
cluster simulations. In this model, each star is give n a certain posi- 
tion and velocity according to the Plummer model jPlummeiligTlh 
with a certain virial radius R\i. The models are assigned virial radii 
of Rm =0.1-1 pc, which are typical for young clusters. Plum- 
mer models are isotropic, and the stellar velocities follow roughly 
a Maxwellian distribution (Fig. [3}. 

(ii) Fractal models. Young sta r cluste r s show a signif- 
icant fraction of substructure (e.g.. iLarsorj 1 199 5: Elmegree: 



2000l: 



Testi et al. 



l2000l : lLada & Ladail2003l : lGutermuth et al,.200a 



eerj 



Allen et al]l2007 ). We set the fractal dimension o- to 1.5 (fractal). 
For comparison, a value a = 3 corresponds to a homogeneous 
sphere with radius R. Each star i s assigned a velocity, as described 
in I Goodwin & WhitworthI ( |2004|) . such that nearby stars have sim- 
ilar velocities. As in the Plummer models, each cluster is assigned 
a radius in the range = 0. 1 - 1 pc. Note, however, the difference 
between the definition of R for the two sets of models. 

The virial ratio Q = -E^lEp of a star cluster is defined as 
the ratio between its kinetic energy Ek and potential energy Ep. 
Clusters with Q = 1/2 are in virial equilibrium, and those with 
2 < 1/2 and Q > 111 are contracting and expanding, respec- 
tively. We study both clusters in virial equilibrium (Q = 1/2), 
as well as clusters with Q = 3/2. The latter value for Q is ex- 
pected f or young clusters with an e ff^ective star forming efficiency 
of 33% jOoodwin & Bastiaij|2006l) . We perform the simulations 
until the clusters are completely dissolved, which is typically of 
the order of 20 - 50 Myr, the timescale at whic h the majority 



of low-mass star clusters ar e destroyed (see, e.g . . iTutukov 



197S; 



2005 



e-g -,L 

iBoutloukos & Lamers|[2003l : iBastian etai]|2005l : [pall et al 

and numerous others). 

3.2 Binarity and multiplicity 

At first, we will consider star clusters that initially consist of sin- 
gle stars only, while later (§ 15. 3t we will also include primordial 
binaries. We do not study the evolution of star clusters with pri- 
mordial higher-order (A' > 3) multiple systems. However, these 
higher-order systems do form in our star cluster simulations. In this 
case, the following three useful quantities describing the multiplic- 
ity of a stellar population can be used: 

S = (B+T + ... )/(5 +B + T + ...) (2) 
= (2S + 3r + ...)/(5 +2S + 3r + ...) (3) 
C = (B + 2T + ... )/(5 +B+T + ...) (4) 

(see, e.g.. lReipurth & ZinneckeJI 19931 : [Kouwenhoven et alj|2005h . 

Here, 5 , B, and T denote the number of single stars, binaries, and 



triples in the system. The quantity S is the multiplicity fraction 
(commonly known as the "binary fraction"). N is the non-single 
star fraction, as 1 - A/ is the fraction of stars that are single. Fi- 
nally, C is the companion star fraction, which describes the average 
number of companions per system, where "system" can refer to a 
single star or multiple system. The number of systems is given by 
S +B + T+..., while 5 -I- 2S -I- 3r -I- . . . denotes the total number of 
individual stars. Clusters that only contain single stars and binary 
systems have S = C. 

3.3 Detection and stability of wide binary and multiple 

systems 

After each simulation, potential binary and multiple systems are 
identified as those pairs with negative energy (see also lParker et al.l 
|2009'). A multiple (A' > 3) system can only survive for a consid- 
erable amount of time if (i) the system is internally stable, and (ii) 
if the outer orbit is stable against perturbations and tidal forces in 
the Galactic field. To ensure internal stability of each level in the 
hierarchy of the multiple system, we impose the Valtonen stability 
criterion Oout/f in > 2st> whe re a^^ and flpm are the semi-major axes 
of the inner and outer orbits. IValtonen et al] J200"8i) find that 



a. «3 1 + 



Mi 



\2/3 

-) (l-e)-"'(^ 



I cos i - cos^ I 



\l/3 

i) , (5) 



where Mi„ is the (total) mass of the inner component, Mout the mass 
of the outer component, ; the relative inclination of the orbits, and e 
the eccentricity of the outer orbit. For a typical system consisting of 
equal-mass stars, a circular outer orbit (e = 0) and a prograde outer 
orbit i = 0, the above expression reduces to 2st ~ 3.7. Systems 
with Oout/fin > 2st are internally stable for at least 10"* revolutions 
of the outer component. For wide binaries, with orbital periods of 
~ 500000 years (flout ~ lO'' au), this corresponds to an internally 
stable period of at least 2 Gyr. 

Wide orbits may additionally be unstable against the tidal 
forces in the Galactic field and interactions with other single stars 
and binaries. We therefore additionally impose a maximum semi- 
major axis of 0.1 pc on the outermost orbit of a binary or multiple, 
motivated by the observed wide binary population (see § The 
stability of wider binaries is difficult to assess. As several binaries 
wider than 0.1 pc are known, our predictions may slightly underes- 
timate the wide binary fraction. The properties of the wide binary 
populations described in this paper therefore pertain to binaries in 
the separation range 10'' au < a < 0.1 pc. Note that these binary 
systems fall well in the category "extremely wide binaries" in the 
IZinneckej jl984) classification of orbital separations. 



4 ANALYTIC AND NUMERICAL ESTIMATES 

Before proceeding to the A^-body simulations in § [5] it is useful 
to first obtain some analytical approximations for the prevalence 
of wide binaries that form during cluster dissolution as well as 
their orbital characteristics. To this end, we first obtain rough esti- 
mates using an analytical approach (§ 14. lb . and subsequently using 
a Monte Carlo approach (§ 14.21 ). 



4.1 Binary formation in Maxwellian velocity space 

Wide binaries may form during the dissociation of a cluster if the 
relative velocity between two stars is sufficiently small that they be- 
come bound once the perturbing cluster potential is removed. Here 
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Figure 3. The distribution of velocities and relative speeds for Plummer 
model with N = 1000, a Kroupa IMF, Q = \I2 and a virial radius of 
i?v = 0.1 pc. Left: The cluster members in the (Vj, i'y)-diagram, where, as 
an example, the encircled pair of stai's represent a potential binary system. 
Right: the distribution of relative speeds (high-resolution histogram) closely 
follows the Maxwellian distribution (solid curve). 




1000 



Figure 4. The number of close neighbours to each star (corresponding to the 
number of (potential) binaries formed per star when A'ncigh * 1) as a func- 
tion of N, v^iit for Plummer models with Ry =0.1 pc. Results are shown 
for models with (2=1/2 (sohd curves) and Q = 3/2 (dashed curves). From 
bottom to top, both sets of curves represent the results for Vcnt = 0.10, 0.15, 
0.20, 0.25, and 0.30 km s"'. 



we assume that a wide binary will form if the relative velocity is 
smaller than or roughly equal to the orbital velocity of a star in 
a wide binary system. For two stars with masses Mi and M2 in a 
circular binary orbit with semi-major axis a, the velocities of the 
individual stars are given by 



and V2 = viq 



(6) 



a(Mi + M2) 

where G is the gravitational constant and q = M2/M1 is the mass 
ratio of the binary system. When adopting, for simplicity, q = I, 
the above expression reduces to: 

\ au/ 



Vorb K 30 km s 



Ml + Ml 



(7) 



where v^i, is the velocity of either of the components. In order to be 
able to form a binary with semi-major axis a, we require velocity 
differences to be smaller than the critical velocity 



For our choice of the IMF, the total mass of a binary system is of 
order 1 Mq. Binaries with a = 3.6 x 10^, 9 x lO'* and 4 x 10"* au 
thus typically require velocity differences of = 0.1, 0.2 and 

0.3 kms"', respectively. 

If the velocity distribution of the stars in a given star cluster 
follows a Maxwell-Boltzmann distribution, then so does also the 
distribution of relative speeds between the stars. We define the rel- 
ative velocity, V = v, - v^, with components (V,-, V,., V;) and mag- 
nitude V = |V|. The distribution over relative speeds is then given 
by: 



1 



teinnev & TremainJ 19871 p. 485), where cr is the one-dimensional 
velocity dispersion. In the Plummer model cr is given by: 



'2=1/2 



(r) = 



I6GMci 



1 + 



I6r 



-1/2 



(10) 



dHeggie & Hu^l2003l) for a cluster in virial equilibrium (i.e. Q = 
I /2). Here, is the total mass of the cluster, and r the distance to 
the cluster centre. We can re-write Eq. dlOt in units more suitable 
for the clusters considered in this paper. First, we set M^i = N{m), 
where A' is the n umber of stars in the cluster and (m) is their average 
mass. Using the iKroupal j200lh IMF (m) = 0.55 Mq. Evaluating 
the velocity dispersion at the (intrinsic) half-mass radius, r = i?i,m ~ 
0.769i?v of the cluster, we find: 



o"e=i/2 = 0.64 x 



0.1 pc 



llOO stars/ 



(11) 



The kinetic energy for a star cluster with 2 = 3/2 is three times that 
of a cluster with Q = 1/2, and therefore the corresponding velocity 
dispersion is 

o"2=3/2 = V3 0-3= 1/2 . (12) 

As an example. Fig. [3] shows the distribution of velocities (v,-, v,.) 
for a Plummer model with A' = 1000 stars, a virial radius R^j = 
0.1 pc, and Q = 1/2. The distribution of relative speeds between 
random pairs of stars in the cluster is shown in the right-hand panel. 
The latter distribution is well approximated by Eq. l[9j with cr = 
1.9 kms"', the velocity dispersion at the half-mass radius is given 
by Eq. 

To find the relative fraction of pairs in a given star clus- 
ter which has a relative speed such that they may become bound 
when the cluster disperses, we integrate Eq. ^ between v = and 
V = Vciit, where v^yit is the critical velocity difference (Eq.[8]l, below 
which we assume that two stars may become bound after cluster 
dissolution: 



£-'P(V)dV /;--exp(-^)K^^V^ 



(13) 



V < 2v„rb = Vr, 



(8) 



where we normalised the fraction to unity by dividing by the inte- 
gral of Eq. ^ between and 00. To find the number of pairs with 
relative speed less than v„it we multiply Fi,hyN- I. Hence: 

/V„eigh = (/V- l)Fi. (14) 

If A/ncish is smaller than unity one might expect that the binary frac- 
tion is proportional to A'neigh- In situations where N„agh is larger (i.e., 
larger than unity) and hence many stars are close to each other in 
velocity space, we might expect to have some competition between 
the stars to stay bound. 

As an example we show in Fig. |4] how N^dO: varies with the 
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Figure 5. Monte-Carlo predictions for tlie dependence of tlie median semi- 
major axis a of tlie escaping wide (10^ au < a < 0.1 pc) binaries {left) 
and wide binary fraction B {right) on the number of stars in a cluster (top) 
and its initial radius R {bottom). Results are shown for Q = i/2 (solid 
curves) and Q = 3/2 (dashed curves). The thin and thick curves correspond 
to the Plummer and fractal models, respectively. The diagonal dotted line 
indicates a = R. The other properties of the modelled clusters are listed in 
TableQ 

number of stars, A^, in clusters with 2=1/2 and Q = 3/2. We show 
results in the velocity range Vcru = 0. 1 - 0.3 kms '. Depending on 
their masses and mass ratios, these velocities correspond to binary 
systems with semi-major axes from ~ 0.1 pc down to ~ 10"* au. 
Velocities of 1.9 kms"' (not shown in Fig.|4j roughly correspond 
to binary systems with a = 10' au. 

The values of A'neigh in Fig. |4] are rather high, as compared to 
wide binary fractions derived in § |4.2| and|5] mainly because /Vneigh 
also contains companion stars outside the separation range 10' au- 
0.1 pc considered throughout this paper. In addition, we believe that 
the predicted values will drop further due to the inefficiency of the 
process. For example, it is not likely (but also not impossible) that 
two stars with nearly the same velocity will form a wide binary 
system, if there are other stars in between them. Furthermore, we 
have only considered the relative velocities, while we have ignored 
the relative positions between the stars. However, this analysis does 
provide a strong upper limit on the (wide) binary fractions we might 
expect after cluster dissolution. 

4.2 Upper limits from a Monte-Carlo approach 

In the previous section we obtained rough estimates for the number 
of newly formed binaries using a Maxwellian velocity distribution. 
However, we were unable to recover the distributions of orbital 
properties, such as the semi-major axis distribution, and we were 
not able to take into account the mass spectrum of stars in the clus- 
ter. We therefore use a somewhat more sophisticated Monte Carlo 
approach, to obtain estimates for these properties as a function of 
cluster size, structure, number of stars, and virial ratio, for the mod- 
els listed in Table [T] Estimates of these properties are obtained us- 
ing an ensemble of initial condition snapshots of each model. 

We identify the potential binaries in each star cluster as fol- 



Figure 6. Same as Fig.|5] but showing the results of the W-body simulations. 



lows. For each star with mass Mi we determine its nearest neigh- 
bour. More precisely, we determine the "most bound" neighbour, 
i.e., the neighbour with mass M2 for which the internal binding en- 
ergy 



1 Ml /Wo , GM1M2 
£h = —r 

2 Ml -I- M2 r 



(15) 



is most negative. Here, r is the distance between the two stars, v 
their velocity difference, and G the gravitational constant. Subse- 
quently, we select those pairs of stars that are each other's mutual 
nearest neighbours, and assume that they will form a binary sys- 
tem with a semi-major axis of approximately r. Note that not all 
of these bound pairs may actually form a binary system, as their 
velocities are perturbed by neighbouring stars. We also ignore the 
possible presence of triple systems and higher-order systems that 
may form. The results, shown in Fig. |5] should therefore be con- 
sidered as a first-order approximation. We will discuss this figure 
in detail in § [5] where we will compare the results with those of 
/V-body simulations (shown in Fig.|6). 

Based on a simple Monte Carlo approach, we find that the 
wide binary fraction decreases with increasing stellar density, and 
mildly decreases with increasing virial ratio. However, several sim- 
plifications have been made, and therefore these results have to be 
interpreted with care. In particular, we have ignored the interaction 
of each star with all other stars; we have ignored two-body interac- 
tions as well as the tidal field of the cluster. In the following section 
we perform a more accurate analysis to obtain the abundance and 
properties of wide binaries formed during cluster dissolution, by 
performing /V-body simulations. We will discuss all properties in 
detail, and compare these to the results obtained using the analyti- 
cal and Monte-Carlo approaches. 



5 RESULTS FROM A'-BODY SIMULATIONS 

The previous two sections have shown that wide binary formation 
during cluster dissolution may well be possible. In particular, small 
dense clusters seem the most likely to form wide binaries. 



The formation of very wide binaries during the star cluster dissolution phase 



Plummer, Q = 1/2 



Plummer, Q = 3/2 



Plummer, Q = 1/2 



Plummer, Q = 3/2 





1.0 




0.8 


g 






0.6 








0.4 


b 


0.2 




0.0 




0.0 0.2 0.4 0.6 0.8 1.0 
Mass ratio 

Fractal, Q = 1/2 




0.0 0.2 0.4 0.6 0.8 1.0 
Mass ratio 



Fractal, Q = 3/2 




0.0 0.2 0.4 0.6 0.8 1.0 
Mass ratio 



0.0 0.2 0.4 0.6 0.8 1.0 
Mass ratio 




0.0 0.2 0.4 0.6 0.8 1.0 
Eccentricity 

Fractal, Q = 1/2 



0.0 0.2 0.4 0.6 0.8 1.0 
Eccentricity 

Fractal, Q = 3/2 




0.0 0.2 0.4 0.6 0.8 1.0 
Eccentricity 



0.0 0.2 0.4 0.6 0.8 1.0 
Eccentricity 



Figure 9. The mass ratio distribution for binary and multiple systems with 
a > 10^ au resulting from the models in Figs.|2]and [8] The solid and dashed 
curves indicate the mass ratio distributions for binaries with Mi > 1.5 Mq 
and M[ < 1.5 Mq, respectively. 



In this section we use A^-body simulations of evolving star 
clusters to study how the properties of the newly formed binaries 
and multiple systems depend on the initial properties of the clus- 
ters. In § 15.11 we describe the orbital properties and multiplicity 
fractions of wide binaries resulting from a typical star cluster. In 
§ l5.2l we show how the results depend on the initial size R and the 
number of stars A' in clusters consisting of initially single stars, for 
the models listed in Table [T] In § l5.3l we study how the results are 
affected by the presence of primordial binaries. 



5.1 Properties of the newly formed binary population 

We perform W-body simulations of Plummer and fractal clusters 
consisting of N = 1000 stars with radii R = OA pc (the virial ra- 
dius for the Plummer models, or the total radius of fractal mod- 
els), with initial virial ratios Q = 1/2 and Q = 3/2 (i.e., models 
P2v, P2e, F2v and F2e in Table[T](. Fifty realisations of each model 
are performed to improve the statistics. The resulting distributions 
over mass, mass ratio, and semi-major axis for the resulting binary 
population are shown in Figs. [7] (Plummer models) and [8] (fractal 
models). 

The left-hand panels in Figs. |7] and [8] show the separation 
distribution f(a) of the resulting binary population. Note that bina- 
ries of all separations are included in these figures, irrespective of 
whether they are actually able to survive in the Galactic field or not. 
The figures illustrate that the separation distribution of the newly 
formed binaries is bimodal, and consists of a small-separation dy- 
namical peak and a large-separation dissolution peak. The tighter 
binaries in the dynamical peak are formed during dynamical en- 
counters in the cluster, and most of them remain mutually bound 
during the further evolution of the cluster. The wide binaries in the 
dissolution peak, on the other hand, are formed during the dissolu- 
tion phase of star clusters. 

The two sets of models with a Plummer density distribution 
result in a small dynamical peak, indicating that dynamical inter- 



Figure 10. The eccentricity distribution f(e) for binary and multiple sys- 
tems with a > lO' au resulting from the models in Figs.|2]and [8] The his- 
tograms indicate f(e), while the curves show the corresponding cumulative 
distributions. The filled circles represent the cumulative thermal eccentric- 
ity distribution f{e) = 2e. 



actions during the lifetime of the cluster generally do not result in 
the formation of close binaries. This is not surprising, as all stars in 
the Plummer models are initially given random velocities. Due to 
the immediate expansion of the Plummer model with Q = 3/2, the 
dynamical peak is completely absent in this case. 

For the models with a fractal density distribution in Fig. [8] 
there are numerous binaries in the dynamical peak. The fact that 
the dynamical peak is stronger for models F2v/F2e than for models 
P2v/P2e is due to both the initial positions and the initial veloc- 
ities being correlated in the fractal models. Although the average 
distance between two random stars is similar for both models, the 
average distance between nearest neighbours in the fractal mod- 
els is smaller (as they are clumpy). As nearby stars in the clumpy 
structure also have similar velocities, frequent dynamical interac- 
tions occur, resulting in a strong dynamical peak. 

For our choice of initial conditions, binaries in the dynami- 
cal peak have separations in the range 1-10^ au, with a median 
value near 50 - 100 au. The median value is set by the typical dis- 
tance between stars in the most densely populated regions of the 
cluster during the formation of these binary systems. Interestingly, 
this also corresponds to the obser ved peak in the Taurus-Auriga bi- 
nary separation distribution (e.g.. lLeinert et alJll993l : [Kroupa et al.l 
1 19991) . 

Binaries in the dissolution peak have a semi-major axis in the 
separation range 10^ au-5 pc. The widest binaries in the dissolution 
peak will immediately break up in the Galactic field, hence our 
choice to study the wide binary population in the separation range 
10^ au - 0. 1 pc throughout this paper. The median separation of 
binaries in the dissolution peak occurs at a « 0.1 - 0.2 pc. As we 
will see later (§ 15.2.2b . this value is set by the initial size of the 
cluster. 

For practical purposes, we consider three ranges in semi- 
major axis: close binaries with a < lO' au, wide binaries with 
10^ au < a < 0.1 pc, and extremely wide binaries with a > 0.1 pc. 
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Figure 7. The semi-major axis distribution (left), the coirelation between mass ratio q and semi-major axis a (middle) and between primary mass and mass 
ratio (right). The histograms in the semi-major axis distribution are normalized such that the maximum value equals unity. The properties of the orbits of 
binary systems and higher-order multiple systems ai'e indicated with the dots and triangles, respectively. For each multiple system with n stellar' components, 
we have included all n - 1 orbits. Results ai'e shown for 50 Plummer models with N = 1000 and R = 0.1 pc, and vrrial ratios of Q = 1/2 (top) and Q = 3/2 
(bottom). The vertical dashed lines indicate a = 10^ au and a = 0.1 pc, respectively. The dashed curve in the right-hand panel indicates the minimum mass 

ratio ^min(Ml) = Mmin/Ml. 



Table 2. The specific binary fraction S for the models shown in Figs. [T] 
and [8] in which the three ranges in semi-major axis ai'e divided with the 
vertical dotted lines. 



Model 


S 


S 


S 


Separation range 


< 10^ au 


U9 au-0.1 pc 


> 0.1 pc 


P2v (N = 1000) 


0.2% 


0.3% 


0.8% 


P2e (N = 1000) 


0.1% 


1.4% 


2.8% 


F2v (N = 1000) 


3.3% 


1.8% 


2.6% 


F2e (N = 1000) 


2.2% 


0.6% 


1.1% 



The limits are indicated with the vertical dotted lines in the figures. 
Most close binaries that are found in star clusters are formed via 
the "normal" star formation process, with the small number seen in 
these simulations formed by dynamical interactions. The wide and 
extremely wide binaries are formed during the cluster dissolution 
phase. Note however, that the vast majority of the extremely wide 
binaries are unstable in the Galactic field, and are ionised quickly 
after their formation. 

For the models in Figs. |7] and [8] the specific binary fraction 
(i.e., the fraction of binary systems in a certain semi-major axis 
range) of the three types of binaries are listed in Table[2] The high- 
est wide binary fractions of a few per cent (in the separation range 
10' au-0.1 pc) are obtained for Plummer models with Q = 3/2, 
and fractal models with Q = 1/2. 



The middle and right-hand panels of Figs. |7] and [8] show 
the correlations between semi-major axis, mass ratio, and primary 
mass, for the binary and multiple (higher-order) systems in each 
of the models. The panels indicate the presence of a large num- 
ber of newly formed multiple systems. These higher-order systems 
are stable in isolation, but a large fraction will not be able to sur- 
vive in the Galactic field, where tidal forces will rapidly remove the 
outer component from the system. Figs. [7] and [8] therefore overes- 
timate the fraction of higher-order multiple systems. Note, in par- 
ticular, the high prevalence of multiple systems in the dynamical 
peak. Many outer components of these systems fall in the dissolu- 
tion peak. These systems are therefore wide higher-order systems. 

For the Plummer models, the dynamical peak consists of sys- 
tems with high masses and high mass ratios. This is a well-known 
signature of mass segregation: the highest-mass st ars sink to the 
cluster centre, where they form close binaries (e.g. 'Heggie & HuJ 
2003). During the dissolution phase of the clusters, these close, 
massive binaries act like single stars when forming a "wide bi- 
nary", which is in fact a wide triple or higher-order multiple system. 
The effect of mass segregation is less visible for the fractal models, 
where dynamical interactions in the subclumps play a greater role. 
However, Fig. [8] still clearly shows that most massive systems are 
mostly close (a < 10^ au) and often higher-order. In addition to 
the triple and higher-order systems formed during the dissolution 
process, s everal higher-order system s may form via dynamical in- 
teractions l lvan den Berk et ani2007l) . 
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Figure 8. Same as Fig. [7] but now for fractal models with N = 1000 and R = 0.1 pc , and virial ratios of Q = 1/2 (top) and Q = 3/2 (bottom); in each case 
fifty realisations have been simulated. 



The correlations between the primary mass and mass ratio dis- 
tributions for the binary and multiple systems are similar to those 
expected from random pairing of individual stars. To first order ap- 
proximation, the masses of the two stars, Mi and M2, in each binary 
are uncorrected; one would therefore expect som ething similar to 
random pairing (e.g., iKouwenhoven et ^|2009|) , where the aver- 
age mass ratio decreases with increasing binary system mass. The 
resulting mass ratio distributions for binary and multiple systems 
with a > 10^ au (i.e., those in the dissolution peak) are shown in 
Fig. 121 which illustrates the dependence of the mass ratio distribu- 
tion on mass. 

Base d on the a nalysis of a sample of 798 common proper mo- 
tion pairs, iTrimbld fl987.) also come to the conclusion that the very 
wide bi nary population i n the field is consistent with random pair- 
ing, andl Valtonei] ( ll997l) come to the same conclusion from their 
simulations of three-body encounters. However, the wide binary 
population does not result from random pairing alone, as the in- 
teraction between two stars depends on their mutual gravitational 
attraction, and the probability of two stars forming a binary is thus 
proportional to the product M[M2. In other words, gravitational fo- 
cusing (e.g., Gaburov et al. 2008) plays an important role. 

Measurements of the eccentricity distribution of wide binaries 
are currently unavailable, due to the large orbital periods and in- 
completeness. If we suspect that the vast majority of wide binaries 
probably have formed dynamically, and as dynamical interactions 
are common among the widest binaries (with respect to closer-in 
binaries), the best guess is pe rhaps the ther mal ecce ntricity distri- 
bution /(e) = 2e (0 < e < 1) ( lHeggielll975l . see also lKroupall2008l 
for a derivation), which results from energy equipartition. The ec- 
centricity distributions resulting for binaries in the dissolution peak 



(a > 10^ au) are shown in Fig.fTO] As expected, the thermal eccen- 
tricity distribution is a good approximation for the newly formed 
binary population. 

If wide binaries form during the dissolution process of a star 
cluster, then the orbital and spin angular momenta of the compo- 
nents should be randomly aligned. On the other hand, if the two 
components formed together in some way it might be that the or- 
bital and spin angular momenta of the components will be corre- 
lated (as s een for example in the observations of ~ 100-au Ae/Be 
binaries bv lBaines et alj2006 l). Therefore, observations of the rela- 
tive alignments of orbital and spin angular momenta could provide 
constraints on the possible formation mechanisms of very wide bi- 
naries. 

Finally, the age diff'erence (between primary and companion 
star) for a populatio n of wide binaries could provide a clue to 
their origin (see, e.g.. lKraus & H illenbrand 2009a). For a star clus- 
ter with a certain age spread, one might expect the components 
of the resulting wide binary population to exhibit a similar age 
diff'erence. This age difference is measurable, but only for young 
(< 10 Myr) binary systems. On the other hand, this age difference 
may be smaller than expected from random pairing, if an initial 
correlation between position and velocity exists. 



5.2 Dependence on cluster properties 

In this section we describe how the properties of the wide binary 
population depend on the initial conditions we assign to a star clus- 
ter, in particular its size R, number of stars A', virial ratio Q, and 
morphology (Plummer sphere or fractal structure). We adopt the 
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cluster properties listed in Table [T] We compare the results that we 
derived earlier using Monte Carlo simulations (Fig.O, with the re- 
sults of A'-body simulations, shown in Fig.|6l 



5.2.1 Dependence on the initial cluster mass 

The top panels of Fig.[6]show the median semi-major axis fln,ed and 
the binary fraction S of wide binaries (10^ au < a < 0. 1 pc) as a 
function of the number of stars W in a cluster. For both the Plummet 
models and the fractal models, a^^^ does not vary significantly with 
A' and the virial ratio Q. The reason for this is that all these models 
have an identical size R. Since R is the most important size scale 
imposed on the modelled star clusters, it determines the size scaling 
(i.e., semi-major axis distribution) of the newly formed binaries. 

The dependence of S on A' and Q is qualitatively the same as 
the analytical predictions shown in Fig.|4]and the Monte Carlo ap- 
proximation shown in Fig.|5]. The wide binary fraction S decreases 
with increasing A' because the stars are further apart in velocity 
space (cf. Fig. O, i.e., the velocity dispersion is larger, and hence 
two neighbouring stars are less likely to form a bound system. 

For the A'-body simulations we find that the fractal model with 
2=1/2 provides the highest wide binary fractions, although the 
difference between models is fairly small (especially when com- 
pared to the difference with increasing A'). Models with Q = 3/2 
generally result in a smaller wide binary fraction than those with 
Q = 1/2, due to the larger distance between the stars in velocity 
space (see Eos. II llandll2t. The curves for the fractal models in 
Figs. [5] and Fig. |6] are almost the same, indicating that the Monte 
Carlo approach provides a good estimate of the wide binary popu- 
lation. For the Plummet models, the Monte Carlo approach predicts 
a binary fraction that is too high, which is due to the fact that the 
positions and velocities of stars in the Plummer models as we ini- 
tialise them are uncorrelated. 



5.2.2 Dependence on the initial cluster size 

The bottom panels of Fig |6] shows the dependence of Onjej and S 
on the initial size R of the clusters. Again, these values are only for 
wide binaries with 10^^ au < a < 0. 1 pc. Note the different defini- 
tions of R: for the Plummer models R represents the virial radius, 
while for the fractal models R represents the radius of the sphere 
enclosing the whole system. Note again the similarity between the 
Monte Carlo approximation shown in Fig.|5]and the A'-body mod- 
els. 

As discussed above, the initial cluster size R determines the 
length scale in each model, and therefore the size scaling of the 
semi-major axis distribution of the newly formed binaries. For ex- 
ample, changing the initial size of the clusters shown in Figs. |7] 
and[8]would simply result in the semi-major axis distribution in the 
left-hand panels being shifted to smaller or larger values of a. 

This direct dependence of f(d) on R is not seen directly in 
Fig. |5] because we only show the results for wide binaries in the 
separation range 10'' au < a < 0.1 pc, and because f{a) is bimodal. 
However, the i?-dependent median semi-major axis and binary frac- 
tion can be explained by the dynamical peak and dissolution peak 
shifting through the range 10^ au < a < 0. 1 pc whilst varying R. 

The highest S is found when either the dynamical peak, or the 
dissolution peak, is centred in the separation range 10^^ au < a < 
0.1 pc. For our choice of the initial conditions, this peak occurs at 
R = 0.025 pc for the Plummer models, when the dissolution peak 
is centred in the range 10^ au - 0. 1 pc. The peak in S occurs atR x 
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Figure 11. The effect of primordial binarity on the formation of wide bi- 
naries, for Plummer models with N = 10 and R = 0.1 pc. The solid and 
dashed curves in each panel indicate the results for Q = 1/2 and Q = 3/2, 
respectively. Top: the effect of a variable primordial binary fraction So- The 
bottom horizontal dotted Une indicates a = 18.2 au, the median semi-major 
axis for primordial binaries. Bottom: the effect of the semi-major axis ag for 
models with a primordial binary frequency So = 50% in which each binary 
has a semi-major axis a = ciq. The dash-dotted lines indicate a = 18.2 au, 
the median semi-major axis of binary systems in the Galactic field. The ver- 
tical dotted line indicates a = lO' au, beyond which all primordial binaries 
are classified as wide binaiies. 



0.15 pc for the fractal models with Q = 1/2 and atR^ 0.6 pc for 
fractal models with Q = 3/2, when the dynamical peak is centred 
in the range 10' au - 0. 1 pc. 

Given our set of initial conditions, compact clusters result in a 
wide binary fraction of 8- 12%, irrespective of virial ratio and mor- 
phology. For more extended clusters, those with a Plummer struc- 
ture and those with a higher virial ratio result in a smaller binary 
fraction. The difference between the Plummer and fractal models 
can be explained by (i) the difference in the definition of R for the 
two sets of models, and (ii) by the different intrinsic separation dis- 
tribution (see the left-hand panels in Figs.|7]and[8j. 

The cluster size R determines the length scale of the system, 
and therefore determines the typical semi-major axis of the newly 
formed wide binaries. Other, less important length scales in the sys- 
tem are the mean distance between two stars, which depends on the 
parameters R, N, and the stellar density distribution (see § 15.2. It . 
as well as the typical semi-major axis of primordial binary systems 
(see §113}. 



5.3 EfTects of primordial binarity 

In the analysis above we have considered star clusters that initially 
consist of single stars only. The results for star clusters with a non- 
zero primordial binary fraction are very similar to the results de- 
scribed above, with the difference that the components of the wide 
"binary" are now in many cases primordial binaries. In other words, 
the majority of the wide "binaries" that formed in the simulations 
described in the previous sections, actually describe the outer orbits 
of wide triple and quadruple systems. 

We predict the properties of these wide triple and quadruple 
systems by performing A'-body simulations of the Plummer models 
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listed in Table[T] but now we include a non-zero primordial binary 
fraction. We perform the simulations with primor dial binary frac- 
tions So ranging from 0% to 100%. We adopt the' Kroupal h995h 
birth period distribution. This distribution is derived from a detailed 
analysis of observed stellar populations, and has the form 



fpiP) = 2.5 ■ 



(logP-logP„„„) 
45 + (log/'-logP„i„)' 



(16) 



for P^,„ < P < Pmax, where logP^,„ = 1, logP^ax = 8.43, and P 
is the period in days. We adopt a thermal eccentricity distribution 
/(e) = 2e(0 < e < 1). We adopt a flat mass ratio distribution f(q) = 
1 with < q = M2/M1 < 1 (i .e., we apply pairing function PCP-I; 
see lKouwenhoven et alj2009b . Subsequently, we generate an initial 
population fro m this bir t h pop ulation, by applying eigenevolution 
as described in iKroupal ( [1993) . All binaries are assigned random 
orientations and orbital phases at the beginning of the simulations. 

Due to the inclusion of binary components, the total mass of 
each cluster increases slightly (up to a maximum of 50% for a pri- 
mordial binary fraction of 100%), although the number of "sys- 
tems", N = S + B remains constant. Strictly speaking, it is thus not 
appropriate to directly compare clusters with and without binaries, 
as we have change d more than one parameter: bin arity and clus- 
ter mass (see, e.g., iKouwenhoven & de Griisll2008h . However, as 
the increase in cluster mass due to adding the companions is rather 
small, we will ignore this issue. 

The results for clusters with a varying primordial binary frac- 
tion So is shown in the top panels of Fig. [TT] For small binary 
fractions, the results are very similar to those of clusters without 
primordial binaries. The properties of the resulting wide binary 
population depend mildly on Q. An increasing Q results in, on av- 
erage, wider binaries, hence in a larger fraction of binaries with 
a > 0. 1 pc, and therefore in a slightly smaller wide binary frac- 
tion. Note that S decreases slightly with increasing Sq. A larger 
primordial binary fraction results in a smaller wide binary fraction, 
possibly because of the destruction of newly formed wide binaries 
by primordial binaries (which have a significantly larger collisional 
cross-section than single stars). 

Whether or not primordial binarity affects the formation of 
wide binaries depends not only on the primordial binary fraction, 
but also on the properties of these binaries: the semi-major axis (or 
period) distribution, the eccentricity distribution and the mass ratio 
distribution. The most important of these is the semi-major axis 
distribution /(a), as it determines the internal binding energy of a 
binary and the cross-section for gravitational interactions between 
binaries and other binaries or single stars. In order to extract the 
dependence on /(a), we simulate clusters in which all binaries have 
a single value for a = ao- We vary ao in each cluster, and determine 
the number of newly-formed binaries. In all simulations we adopt a 
primordial binary fraction of 50%, a flat mass ratio distribution and 
a thermal eccentricity distribution. 

The results of these simulations are shown in the bottom pan- 
els of Fig.[TT] For models with ao < 10^ au, most primordial bina- 
ries survive, while additional wide binary, triple, and quadruple sys- 
tems are formed. In fact, the resulting wide binary fraction is practi- 
cally independent of the primordial binary fraction Sq. For models 
with ao > 10^ au, all primordial binary systems are classified as 
wide binaries. For these models we therefore have S ~ So and a 
median semi-major axis equal to ao, which results in the glitches at 
a = 10^ au in Fig.fTTI 

The wide orbits are part of systems with 2, 3, and 4 compo- 
nents. They are formed by randomly pairing single stars and pri- 
mordial binary systems together. The number of multiple systems 



of each degree can thus be estimated by simply calculating the 
probability of randomly drawing a single-single, single-binary, and 
binary-binary pair. When assuming a primordial binary fraction So, 
the multiplicity distribution of the resulting wide population can be 
estimated as follows: 



Wide binary fraction : 
Wide triple fraction : 
Wide quadruple fraction : 



S(l-So)' 
2SSo(l-So) 
SSL 



(17) 
(18) 
(19) 



where we have made the assumption that none of the primordial 
binary systems has broken up. 

All models shown in Fig. [TT] result in wide binary fractions 
S ~ 8% that are more or less independent of So and ao. The value 
of S is therefore primarily determined by the initial values of the 
number of system N in the cluster, and its initial size R. 

If we assume that the wide orbits in the bottom panels of 
Fig. [TT] (where So = 50%) are formed of randomly paired com- 
ponents (i.e., single stars or primordial binaries), we can calculate 
the fraction of higher-order multiple systems among the S = 8% 
wide binaries. Among these, we predict that 25%, 50%, and 25%, 
are binary, triple, and quadruple systems, respectively. In this exam- 
ple, we thus expect 75% of the "wide binaries" to be higher-order 
multiple systems. Due to the random process, the outer orbits of 
these systems are expected to be uncorrelated with the inner orbits 
or stellar spin axes. 

The ratios between wide binary, triple, and quadruple systems 
are therefore indicative of So. A survey for higher-order multiplic- 
ity among "wide binary systems" can thus be used to constrain 
the primordial binary fraction. Given the fact that the majority of 
stars do form in binary systems, we predict a very high fraction of 
higher-order multiple systems among wide "binary" systems; see 
Fig. [12] Our proposed mechanism could explain the exi stence of 
the observed wide multiple systems iMam aiek et alJ2009l), and our 
predictions are strongly supp orted by the surveys of Makarov et al.l 
( 2008) and Fahertv et al] HoiOl) . who find that a high fraction of the 
common proper motion pairs in their survey contain inner binary or 
triple systems, which is significantly higher than in populations of 
other types of binaries. 



6 SUMMARY AND DISCUSSION 

Observations have shown that ~ 15% of binaries are wide (a > 
10' au). These wide binaries are difficult to explain as being the 
result of star formation as it is difficult to see how wide binaries 
can form (especially those > 10* au), and they would be rapidly 
destroyed in clustered star forming environments. Whilst 10 - 30% 
of stars do appear to form in an 'isolated' envirormient in which 
such binaries could possibly survive, in order to explain the fraction 
of wide binaries, almost all stars forming in isolated environments 
would have to form wide binaries. Further, such wide binaries can- 
not be formed later in any significant numbers by dynamical inter- 
actions in the Galactic field. 

In this paper we study the possibility of wide binary formation 
during the dissolution phase of star clusters, in particular, during the 
rapid expansion of clusters after gas expulsion. We study this pos- 
sibility using (1) an analytical approach in an idealised situation, 
(2) a Monte Carlo approach, and (3) detailed A'-body simulations. 
Our main conclusions are as follows: 

(i) The wide binary fraction S among the dissolved stellar popu- 
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Figure 12. Given the fact that most star form in binary systems, it is ex- 
pected that the majority of the wide "binary" systems are in fact part of a 
triple or quadruple system. 



lation ranges between 1% and 30%, depending on the cluster prop- 
erties. 

(ii) More massive star clusters result in a smaller wide S than 
low-mass clusters. Clusters with a spherical, smooth stellar density 
distribution form fewer wide binaries than substructured clusters of 
the same size and mass. This is due to the fact that the average dis- 
tance between nearest neighbours is smaller for substructured clus- 
ters. Expanding (post-gas expulsion) star clusters produce a larger 
S than those starting out of equilibrium. 

(iii) The typical semi-major axis a of the newly formed bina- 
ries is similar to the initial size R of the star cluster from which 
they were born. The resulting semi-major axis distribution is gen- 
erally bimodal, consisting of a dynamical peak with binary systems 
formed by dynamical interactions, and a dissolution peak with bi- 
nary systems formed during the cluster dissolution phase. 

(iv) The formation of wide binaries during the star cluster disso- 
lution phase is a random process, resulting in the following orbital 
properties. The eccentricity distribution of the wide binaries is ap- 
proximately thermal: /(e) a 2e for < e < 1. The mass ratio dis- 
tribution of the wide binaries is the result of gravitationally-focused 
random pairing. In a wide binary, the orbital and spin angular mo- 
menta are uncorrelated. 

(v) Star clusters with a non-zero primordial binary population 
form wide triple and quadruple systems, i.e., the components of a 
newly-formed wide "binary" can themselves be close primordial 
binaries, rather than single stars. The ratio of triple to quadruple 
systems among very wide orbits is therefore indicative of the pri- 
mordial binary fraction Sq. Given that So is large, we predict a 
high frequency of triple and quadruple systems among the known 
wide "binary" systems, which is supported by existing surveys for 
higher-order multiplicity among wide binary systems. 

Throughout this paper we have made predictions of the prop- 
erties of the wide binary population resulting from the dissolution 
of individual clusters. In order to compare our results with obser- 
vations, we should therefore take into account the fact that the field 
star population is made up of the stars resulting from an ensemble 
of clusters of different sizes and masses. The initial cluster mass 



distr ibution may be approximated by /(M) oc M ^ with y a 2 (see , 
e.g..lZhang &FaI]||l999l : lAshman & ZepjbOOll : iBik etalj[2003l : 
iHunter et alj l2003h . Given the number of stars A' = M^Jim), 
where (m) is the aver age mass of a sta r, this distribution is equiv- 
alent to f{N) oc A^-''. lOev etai] j2004l) suggest that the above ex- 
pression can be extrapolated down to A^mm = 1- The upper limit 
for the initial cluster mass distribution is M^ax ~ 10^ M© (e.g., 
Ide Griis & Parmentiei^ '2007, and references therein). The resulting 
binary fraction Bf for the ensemble of stars (i.e., the field star pop- 
ulation) is then given by: 



(20) 



where S(A'd) is the cluster mass dependent wide binary fraction. 
The numerator in the above expression is proportional to the num- 
ber of binaries, and the denominator is proportional to the to- 
tal number of stars in the ensemble of clusters. In addition, the 
size and dissolution time of a star cluster, and therefore the wide 
binary fraction, may also depend on its Galactic location (e.g., 
^Baumgardt & Makino| l2003). An inspection of Fig.|6]shows that an 
extrapolation toN x 10* results in a wide binary fraction of several 
per cent; smaller than the observed 15%, irrespective of the choices 
for R, Q, and the morphology of the cluster. Although we predict 
rather small values, our back-of-the-envelope calculation does re- 
sult in the right order of magnitude for the wide binary fraction in 
the Galactic field. It is clear, however, that a deeper investigation 
is required to accurately recover the properties of the wide binary 
population in the field. In particular, a wider range of star cluster 
morphologies has to be considered, by varying the fractal dimen- 
sion and position-velocity correlations of individual star clusters. 

Our proposed formation mechanism for very wide binaries 
predicts at least several common proper motion pairs in and around 
dissolving star clusters and moving groups. For example, the mech- 
anism may well explain the presence of the th ree co mmon prope r 
motion pairs in the moving groups studied bv lciarke et al] ( l2009h . 
The future prospects in wide binary research are bright: an enor- 
mous number of wi de binaries are expected to be found with 
the GAIA missioiH (Perrvman et al." 200ll ; lTuron et alj|200 5'l and 
LAMOST0 (Chu 1998; Stone 2008). These datasets should help 
determine the true fraction of wide binaries and their orbital pa- 
rameters. 



ACKNOWLEDGEMENTS 

We would like to thank Anthony Whitworth (our referee) and Si- 
mon Portegies Zwart for a useful discussion on this topic. MBNK 
was supported by the Peter and Patricia Gruber Foundation through 
the PPGF fellowship, the Peking University One Hundred Tal- 
ent Fund (985), and by PPARC/STFC (grant PP/D002036/1). We 
would like to acknowledge CICS for the provision of research 
computing facilities through the White Rose Grid. We acknowl- 
edge research support from and hospitality at the International 
Space Science Institute in Bern (Switzerland), as part of an Inter- 
national Team programme. The authors acknowledge the Sheffield- 
Bonn Royal Society International Joint Project grant, which pro- 
vided financial support and the collaborative opportunities for this 



http://www.esa.inl/science/gaia 
http://www.lamost.org 



The formation of very wide binaries during the star cluster dissolution phase 13 



work. The calculations performed by DM were carried out on com- 
puter hardware purchased from the Royal Physiographic Society in 
Lund. RJP acknowledges financial support from STFC. 



REFERENCES 

Allen C, Poveda A., Hemdndez-Alcantara A., 2007, in Hartkopf 
W. I., Guinan E. R, Harmanec P., eds, lAU Symposium Vol. 240 
of lAU Symposium, Halo Wide Binaries and Moving Clusters as 
Probes of the Dynamical and Merger History of our Galaxy, pp 
405^13 

Allen C, Poveda A., Herrera M. A., 2000, A&A, 356, 529 
Allen L., Megeath S. T., Gutermuth R., Myers P C, Wolk S., 
Adams F. C, MuzeroUe J., Young E., Pipher J. L., 2007, in 
Reipurth B., Jewitt D., Keil K., eds, Protostars and Planets V The 
Structure and Evolution of Young Stellar Clusters, pp 361-376 
Ashman K. M., Zepf S. E., 2001, AJ, 122, 1888 
Bahcall J. N., Hut P, Tremaine S., 1985, ApJ, 290, 15 
Bahcall J. N., Soneira R. M., 1981, ApJ, 246, 122 
Baines D., Oudmaijer R. D., Porter J. M., Pozzo M., 2006, MN- 
RAS, 367, 737 

Bastian N., Gieles M., Lamers H. J. G. L. M., Scheepmaker R. A., 

de Grijs R., 2005, A&A, 431, 905 
Baumgardt H., Makino J., 2003, MNRAS, 340, 227 
Bik A., Lamers H. J. G. L. M., Bastian N., Panagia N., Romaniello 

M., 2003, A&A, 397, 473 
Binney J., Tremaine S., 1987, Galactic dynamics. Princeton, NJ, 

Princeton University Press, 1987 
Boutloukos S. G., Lamers H. J. G. L. M., 2003, MNRAS, 338, 

717 

Burgasser A. J., Kirkpatrick J. D., Reid 1. N., Brown M. E., 

Miskey C. L., Gizis J. E., 2003, ApJ, 586, 512 
Caballero J. A., 2009, A&A, 507, 251 
Chaname J., Gould A., 2004, ApJ, 601, 289 
Chandrasekhar S., 1944, ApJ, 99, 54 
Chu Y, 1998, Highlights of Astronomy, 11, 493 
Clarke J. R. A., Pinfield D. J., Galvez-Ortiz M. C, Jenkins J. S., 

Bumingham B., Deacon N. R., Jones H. R. A., Pokomy R. S., 

Barnes J. R., Day- Jones A. C, 2009, MNRAS, p. 1923 
Close L. M., Richer H. B., Crabtree D. R., 1990, AJ, 100, 1968 
de Grijs R., Parmentier G., 2007, Chinese Journal of Astronomy 

find Astrophysics, 7, 155 
Duquennoy A., Mayor M., 1991, A&A, 248, 485 
Eggers S., Keller H. U., Kroupa P., Markiewicz W. J., 1997, 

Planet. Space Sci., 45, 1099 
Ehnegreen B. G., 2000, ApJ, 530, 277 
ESA 1997, VizieR Online Data Catalog, 1239 
Faherty J. K., Burgasser A. J., West A. A., Bochanski J. J., Cruz 

K. L., Shara M. M., Walter F M., 2010, AJ, 139, 176 
Fall S. M., Chandar R., Whitmore B. C, 2005, ApJ, 631, L133 
Fischer D. A., Marcy G. W., 1992, ApJ, 396, 178 
Gaburov E., Gualandris A., Portegies Zwart S., 2008, MNRAS, 

384, 376 

Garnavich P M., 1988, ApJ, 335, L47 
Goodman J., Hut R, 1993, ApJ, 403, 271 
Goodwin S. R, 2009, ArXiv:091 1.0795 
Goodwin S. P, Bastian N., 2006, MNRAS, 373, 752 
Goodwin S. R, Kroupa P, 2005, A&A, 439, 565 
Goodwin S. P., Kroupa P., Goodman A., Burkert A., 2007, in 
Reipurth B., Jewitt D., Keil K., eds, Protostars and Planets V 



The Fragmentation of Cores and the Initial Binary Population, 
pp 133-147 

Goodwin S. R, Whitworth A. R, 2004, A&A, 413, 929 
Gould A., Bahcall J. N., Maoz D., Yanny B., 1995, ApJ, 441, 200 
Gould A., Eastman J., 2006, astro-ph/06 10799 
Gutermuth R. A., Megeath S. T, Pipher J. L., Wilhams J. P., Allen 
L. E., Myers R C, Raines S. N., 2005, ApJ, 632, 397 
Hartigan R, Strom K. M., Strom S. E., 1994, ApJ, 427, 961 
Heggie D., Hut P., 2003, The Gravitational Million-Body Prob- 
lem: A Multidisciplinary Approach to Star Cluster Dynamics. 
Cambridge University Press, 2003 
Heggie D. C, 1975, MNRAS, 173, 729 
Hernandez X., Lee W. H., 2008, MNRAS, 387, 1727 
Hills J. G., 1975, AJ, 80, 809 

Hunter D. A., Elmegreen B. G., Dupuy T. J., Mortonson M., 2003, 

AJ, 126, 1836 
Jiang Y, Tremaine S., 2009, MNRAS, p. 1640 
Kobulnicky H. A., Fryer C. L., 2007, ApJ, 670, 747 
Kouwenhoven M. B. N., Brown A. G. A., Goodwin S. P., Porte- 
gies Zwart S. F, Kaper L., 2009, A&A, 493, 979 
Kouwenhoven M. B. N., Brown A. G. A., Portegies Zwart S. R, 

Kaper L., 2007, A&A, 474, 77 
Kouwenhoven M. B. N., Brown A. G. A., Zinnecker H., Kaper L., 
Portegies Zwart S. F, 2005, A&A, 430, 137 
Kouwenhoven M. B. N., de Grijs R., 2008, A&A, 480, 103 
Kraus A. L., Hillenbrand L. A., 2009a, ApJ, 704, 531 
Kraus A. L., Hillenbrand L. A., 2009b, ApJ, 703, 1511 
Kroupa P, 1995, MNRAS, 277, 1507 
Kroupa R, 2001, MNRAS, 322, 231 

Kroupa R, 2008, in S. J. Aarseth, C. A. Tout, & R. A. Mardling 

ed., Lecture Notes in Physics, Berlin Springer Verlag Vol. 760 of 
Lecture Notes in Physics, Berlin Springer Verlag, Initial Condi- 
tions for Star Clusters, p. 181 

Kroupa R, Burkert A., 2001, ApJ, 555, 945 

Kroupa P., Petr M. G., McCaughrean M. J., 1999, New Astron- 
omy, 4, 495 

Lada C. J., 2006, ApJ, 640, L63 

Lada C. J., Lada E. A., 2003, ARA&A, 41, 57 

Larson R. B., 1995, MNRAS, 272, 213 

Latham D. W., Schechter P., Tonry J., Bahcall J. N., Soneira R. M., 

1984, ApJ, 281, L41 
Leinert C, Zinnecker H., Weitzel N., Christou J., Ridgway S. T, 

Jameson R., Haas M., Lenzen R., 1993, A&A, 278, 129 
Lepine S., Bongiomo B., 2007, AJ, 133, 889 
Longhitano M., Binggeh B., 2010, ArXiv:1001.3819 
Makarov V. V, Zacharias N., Hennessy G. S., 2008, ApJ, 687, 566 
Mallada E., Fernandez J. A., 2001, in Revista Mexicana de As- 

tronomia y Astrofisica Conference Series Vol. 11 of RMAA 

Conference Series, Dynamical Evolution of Wide Binaries, p. 27 
Mamajek E. E., Kenworthy M. A., Hinz P. M., Meyer M. R., 2009, 

ArXiv:09 11.5028 
Mason B. D., Gies D. R., Hartkopf W. I., Bagnuolo Jr. W. G., ten 

Brunmielaar T., McAlister H. A., 1998, AJ, 115, 821 
Mengel S., Lehnert M. D., Thatte N., Genzel R., 2005, A&A, 443, 

41 

Oey M. S., King N. L., Parker J. W., 2004, AJ, 127, 1632 

Opik E., 1924, Tartu Obs. Publ., 25, No. 6 

Parker R. J., Goodwin S. P., Kroupa P., Kouwenhoven M. B. N., 

2009, MNRAS, 397, 1577 
Ferryman M. A. C, de Boer K. S., Gilmore G., H0g E., Lattanzi 

M. G., Lindegren L., Luri X., Mignard R, Pace O., de Zeeuw 

R T., 2001, A&A, 369, 339 



14 M. B. N. Kouwenhoven, S. P. Goodwin, R. J. Parker, M. B. Davies, D. Malmberg and P. Kroupa 



PlummerH. C, 191 1, MNRAS, 71, 460 

Portegies Zwart S. R, McMillan S. L. W., Hut P., Makino J., 2001, 
MNRAS, 321, 199 

Poveda A., Allen C, 2004, in Allen C, Scarfe C, eds. Rev. Mexi- 
cana Astron. Astrofis. Conference Scries The distribution of sep- 
arations of wide binaries of different ages. p. 49 

Poveda A., AUen C, Hernandez-Alcantara A., 2007, in lAU Sym- 
posium Vol. 240 of LAU Symposium, The Frequency Distribu- 
tion of Semimajor Axes of Wide Binaries: Cosmogony and Dy- 
namical Evolution, pp 417^25 

Quinn D. P., Wilkinson M. I., Irwin M. J., Marshall J., Koch A., 
Belokurov V., 2009, MNRAS, 396, Lll 

Radigan J., Lafreniere D., Jayawardhana R., Doyon R., 2009, ApJ, 
698, 405 

Reipurth B., Zinnecker H., 1993, A&A, 278, 81 

Retterer J. M., King I. R., 1982, ApJ, 254, 214 

Scally A., Clarke C, McCaughrean M. J., 1999, MNRAS, 306, 

253 

Scholz R.-D., Kharchenko N. V., Lodieu N., McCaughrean M. J., 

2008, A&A, 487, 595 
Shatsky N., Tokovinin A., 2002, A&A, 382, 92 
Soderhjelm S., 2007, A&A, 463, 683 
Stone R., 2008, Science, 320, 34 

Testi L., Sargent A. 1., Olmi L., Onello J. S., 2000, ApJ, 540, L53 

Trimble V., 1987, Astronomische Nachrichten, 308, 343 

Turon C, O'Flaherty K. S., Perryman M. A. C, eds, 2005, The 

Three-Dimensional Universe with Gaia Vol. 576 of ES A Special 

Publication 
Tutukov A. v., 1978, A&A, 70, 57 

Valtonen M., Myllari A., Orlov V., Rubinov A., 2008, in lAU 
Symposium Vol. 246 of lAU Symposium, The Problem of Three 
Stars: Stability Limit, pp 209-217 

Valtonen M. J., 1997, ApJ, 485, 785 

van Albada T. S., 1968, Bull. Astron. Inst. Netherlands, 20, 47 
van den Berk J., Portegies Zwart S. R, McMillan S. L. W, 2007, 

MNRAS, 379, 111 
Vereshchagin S., Tutukov A., lungelson L., Kraicheva Z., Popova 

E., 1988, Ap&SS, 142, 245 
Ward-Thompson D., Andre P., Crutcher R., Johnstone D., Onishi 
T., Wilson C, 2007, in B. Reipurth, D. Jewitt, & K. Keil ed., 
Protostars and Planets V An Observational Perspective of Low- 
Mass Dense Cores 11: Evolution Toward the Initial Mass Func- 
tion, pp 33^6 
Wasserman I., Weinberg M. D., 1991, ApJ, 382, 149 
Weinberg M. D., Shapiro S. L., Wasserman L, 1987, ApJ, 312, 
367 

Zhang Q., Fall S. M., 1999, ApJ, 527, L81 

Zinnecker H., 1984, Ap&SS, 99, 41 

Zinnecker H., Yorke H. W, 2007, ARA&A, 45, 481 

This paper has been typeset from a T^(/ LJTb^ file prepared by the 
author. 



